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ABSTRACT 

We present time-dependent simulations of a two-phase accretion flow around a black hole. 
The accretion flow initially is composed of an optically thick and cool disc close to the mid- 
plane, while on top and below the disc there is a hot and optically thin corona. We consider 
several interaction mechanisms as heating of the disc by the corona and Compton cooling of 
the corona by the soft photons of the disc. Mass and energy can be exchanged between the 
disc and the corona due to thermal conduction. For the course of this more exploratory work, 
we limit ourselves to one particular model for a stellar mass black hole accreting at a low 
accretion rate. We confirm earlier both theoretical and observational results which show that 
at low accretion rates the disc close to the black hole cannot survive and is evaporated. Given 
the framework of this model, we now can follow through this phase of disc evaporation time 
dependently. 

Key words: black hole physics - galaxies: jets - X-Rays: binaries. 



1 INTRODUCTION 



X-Ray spectral observations of accretion powered black holes in bi- 
nary systems show that there are at least two spectral components 
present. The lower energy, thermal component is usually attributed 
to the geometrically thin but optically thick disc, which emits 
as a black-body or modified black-body (for the latest spectra l 
model, BHSPEC, see iDavi s et aP l2005l : loavis & Hubenvl 120061) . 
The higher energy, power-law component reflects the presence of 
a possibly geometrically thick but optically thin, hot corona. The 
relative importance of these components changes with time, espe- 
cially during periods of strong variability when these sources un- 
dergo outbursts. If the power-law component dominates, the state 
is called low/hard (LH) state, while if the black-body dominates the 
state is called high/soft (HS) state. The nomenclature 'hard/soft' re- 
flects the shape of the spectrum, while 'high/low' reflects the pre- 
su med accretion rate in these states. For a review on spectral states, 
see iMcClintock & RemillardI (l2004 . 

In general this picture holds for high-mass black holes, i.e. 
AGN (active galactic nuclei), as well. However the low-mass black 
holes, i.e. X-Ray Binaries, are the best studied objects with regards 
to the temporal behaviour There the timescale for outbursts is on 
the order of weeks to months and years, while in AGN it is much 
longer. Hence theoretical models of accretion flows around black 
holes concerning the long-term stability properties are observation- 
ally best tested against X-Ray Binaries. 

Many authors have addressed the thermal and physical states 
of the accretion flow in these systems. Much of the work has con- 
centrated on what is required of the thermal and physical states of 



the accretion flow in order to reproduce the X-ray spectra. In these 
papers, the system is usually assumed to be in a steady configura- 
tion, and the effort is devoted to computation of the resulting spec- 
trum. There has also been consideration given to physical and ther- 
mal interactions between the proposed cool (disc) and hot (corona) 
plasmas, and it is these considerations which have the strongest 
bearing on the time-dependent behaviour of the flow. We briefly 
discuss some of this work. 

Eardley et al. ( 1975) and IShapiro et al.l (1 19761) introduced the 
concept of a two-phase accretion model. Their disc consists of an 
inner, hot, optically thin and geometrically thick corona, and an 
outer, standard, optically thick, but geometrically thin, disc. The 
plasma in the inner parts of the corona is a two- temperature plasma, 
where electrons and ions interact via Coulomb collisions. The main 
aim of the model was to provide an explanation of the spectrum of 
CygX^L 

IWandel & Um] (Il99lh and IWandel & Liatigl (Il99lh present 
models for spectra of Seyfert galaxies given three different underly- 
ing assumptions on the flow geometry. They discuss a "sandwich" 
geometry where the corona and the disc literally for m a sandwich. 
Coron a and disc are either separated radially as in 'Shapiro et alj 
(1 19761) or vertically. Models similar to Wandel & Liang (1991) 
were developed by Haardt & Maraschi| ( 1 1 99 ll. 1 19931) . 

While these models were mostly concerned with the expla- 
nation of the observed spectra of the corresponding object, others 
considered the structural properties of these two-phase discs and 
the interactions between the phases. 

IZvcki et all ( 1 1995b and I Witt et all ( 1 19971) present a two-phase 
accretion disc model where the two phases interact by Compton 
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cooling of the corona by soft photons and the disc and heating of 
the disc by illumination of the corona, respectively. They carry out 
a extensive parameter study of their model, and also include the 
effe cts of mass and energy loss t o a wind. 

iBegelman & McKeg ( Il990l) examine effects of thermal con- 
duction on two-phase media. They consider the thermal interaction 
of the hot and cold phases under isobaric conditions and develop 
criteria to determine which of the the hot or cold phase shrinks or 
growS; 

iRozariska & Czernvl jlOOOal) discuss radiative and conductive 
equilibrium for two-phase accretion disc models. They find stabil- 
ity/instability strongly depen dent on the assumed coron al heating 
mechanism. In a later paper ( Rozanska & Czernv'2000b') they con- 
sider mass loss/gain of both phases due to thermal conduction for 
stationary accretion disc m odels following the considerations of 
IBegelman & McKee" ('l99tf). 

Dullemond ( 1999) presents an semi-analytic model of disc 
evapo ration by therma l conductio n. His model is based on 
iMever & Mever-Hofmeisterl ( Il994l) . He calculates the detailed 
structure of the transition layer between disc and corona for an 
one-temperature corona. In order to dissolve the disc, he needs a 
large evaporation rate which only can be achieved if the conductive 
length scale is larger than the radius, so that radial thermal con- 
duction enters the problem. He gives caveats and questions to what 
extent this type of accretion flow can be represented in a ID model. 

One important point in all two-phase models discussed above 
is how to determine the relative amounts of luminosity produced in 
each of the two phases. Most authors simply introduce a parameter 
which predetermines the fraction of lu minosity produced or or the 
fraction of accretion rate in each phase. iNakamura &Osakil ( fl993h 
assume viscous heating to operate in both the hot corona and cool 
disc. Given their constraints on the scaleheight of both parts, they 
find that most of the heat is produced in the corona. Some (e.g. 
IZvcki et al.ll995h constrain this important parameter through some 
boundary conditions at the interface between corona and disc. All 
of these models however assume a stationary accretion flow. 

It was realised fairly early that accretion in a standard 
optically thick disc around a black hole is viscously unstabl e 
at near-Eddington accretion rates dLightman & EardlevI Il974l) . 
Then electron sc a ttering renders the accretion flow unstable. 
lAbramowicz et al] (Il988l) discovered that in these circumstances 
the radial advection of energy cannot be neglected. They dubbed 
such discs "slim discs" and found that these instabilities are sta- 
bilised by radial advection of energy. The name comes from the fact 
that the ratio of disc scale-height to radial distance, H/R, is close 
to unity (thermal energy close to kinetic energy), so that the usual 
'thin disc' approximation is at best marginally applicable. These 
accretion flows are optically thick ADAFs (advection dominated 
accretion flow). These flows are applicable to the high/soft state. 

A simila r pictu re h olds for optically thin accretion flows. 
iNaravan & Yil ( 1 1995b and llchimarul ( 1 19771) proposed an optically 
thin counterpart to the optically thick ADAF. They consider an op- 
tically thin disc, i.e. corona, close to the black hole and an standard 
optically thick disc further out. Their model is applicable to lower 
accretion rates or the so-called low/hard state. The exact transition 
rate from the optically thin ADAF to the optically thick disc is not 
predicted by these models, but is determined to match the observa- 
tions. 

Many of the observed black hole accretion discs show state 
transitions, i.e. tr ansitions from the high/soft to th e low/hard state 
and vice versa jMcClintock & RemillardI l2004l) . The spectrum 
changes from very soft and black-body dominated (attributed to 



the disc) to a hard, power law dominated spectrum (attributed to 
the corona). These changes appear to reflec t changes in the a ccre- 
tion rate and flow geometry as discussed bv lEsin et alj (Il997h . 

In this paper we introduce the various physical concepts which 
are required in order to determine the time-dependence of such ac- 
cretion flows. Many of the physical properties of such flows are 
not well-defined and it is necessary to make many simplifying as- 
sumptions along the way. We shall need to consider three possible 
disc states, when the disc is a standard disc with little or no corona, 
when the disc has a sizeable corona, and when the disc evaporates 
altogether and there is only a corona present. We also need to de- 
termine how the accretion flow transitions from one of these states 
to another. For example, if the disc is truncated at low accretion 
rates and the inner flow is filled with a hot corona, then we need to 
determine the transition radius below which the disc cannot exist 
any more. Thus we need to introduce a time-dependent two-phase 
model for black hole accretion discs, where mass can be exchanged 
between the two phases owing to thermal conduction. 

Because the modelling of all these processes is computa- 
tionally demanding (apart from all the physical complications 
this comes about basically because we need to follow thermal 
timescales in the inner hot disc for several viscous timescales of 
the outer cool disc), we present here just one example of a time- 
dependent computation. We take as an initial condition a steady- 
state flow where the hot and optically thick corona sandwiches a 
cool and geometrically thin disc for which the accretion rates in 
each component are predetermined. This is similar to some of the 
models derived in the literature of the kind invoked to explain the 
spectral characteristics of these objects (see above). We then fix 
the total accretion rate at the outer boundary and allow the disc 
to evolve, permitting radial viscous evolution of both the disc and 
coronal components, and well as thermal energy and mass flow be- 
tween them. 

We give a detailed description of the physical assumptions we 
make to describe model in Sect.|2] These assumptions of necessity 
over-simplify the problem and are certainly not unique, they should 
however give a reasonable indication of the sort of behaviour we 
might expect from such discs. In Sect. [3] we discuss the relative 
contributions of the corona and disc to the total luminosity and we 
also briefly review some observational results. We then discuss the 
time-dependent equations for the two-phase flow in Section |4] and 
describe the numerical setup in Sect. [5] We describe our results in 
Sect. [6] discuss them in Section|7]and end with the conclusions in 
Sect. [8] 



2 THE MODEL 

We assume that at a given radius the general accretion flow consists 
of an optically thick but geometrically thin accretion disc which is 
sandwiched by a hot geometrically thick and optically thin corona 
above and below. In the following we refer to these parts of the flow 
as disc and corona, respectively. Note that at any radius, it may be 
that either the disc, or the corona, is not present. 

We describe the physics of the disc in Section lTTI the corona 
in Sect. 12.21 and the transition layer in 12.31 The disc or the corona 
might cease to exist if either of them is evaporated, we discuss the 
treatment of these cases in Section |2!4l 
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2.1 The disc 

We model the disc as a standard accretion disc in the so-called one- 
zone approximation. Hence the density pd and temperature just 
represent some average values at a given radius R. Taking constant 
density, the equation of hydrostatic equilibrium gives a pressure 
distribution of the form 
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(1) 



where is the rotational frequency of the disc which we assume to 
be Keplerian, i.e. = \/ GM/R^. Pd{z) is the pressure at a height 
z above the midplane. The pressure in the midplane, which we also 
take to be the representative pressure for the disc, contains contri- 
butions of gas and radiation pressure, 



/'d = /'d(0) = Pd— + 



3c 



d ' 



(2) 



where we consider a fully ionised, solar metallicity gas. Thus we 
set the mean molecular weight ^ = 0.6. Here k is the Boltzmann 
constant and a is the Stefan-Boltzmann constant. 

When the disc is optically thick, the disc radiates as a black 
body and we can treat the radiative transfer in the diffusion approx- 
imation, i.e. we write for the one-sided radiative losses in the disc 

4a A 
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(3) 
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where Kr is the Rosseland mean of the op acity. We use the 
tabulated values of the OPAL opacity proiectfeogers & Ig lesiasI 
Il992; Iglesias & Rogers' 1 19961) for the solar composition of 
iGrev esse & Noels ( 1993). We use the X = 0.7 set of their opac- 
ity tables {X is the hydrogen mass fraction of the matter) for solar 
metallicity (Z = 0.02). 

Energy in the disc is generated by viscosity v. We use the 
parametrisation bv .Shaku ra & Sunvaev ( 1973.) . i.e. 



ac.H 



(4) 



where a is the viscosity parameter (typically smaller than unity, we 
use a = 0. 1), Cs ~ \/P /p the sound speed in the medium and H the 
scale-height. The viscosity produces a torque 
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which gives rise to viscous dissipation 
9 



Qt 



rVdPd^d^^d • 



(5) 



(6) 



In both cases we made use of the fact that the disc rotates at the 
Keplerian speed, i.e. 3(lognd)/9(log7?) = —3/2. 

2.2 The corona 

As for the disc, we take again density within the corona constant 
and equal to some representative value. The pressure distribution 
in the corona at height x = above the disc/corona interface 

is given by 



Pc{x + Hd) = Pc{Hd) - pc^t [Hd + ^x]x 



(7) 



where Pc is the (representative) density in the corona, x the height 
above the disc and the height of the disc. The point where Pc {x+ 
Hd) = defines the thickness of the corona, i.e. 



1 



//c=0, 



(8) 



We set Pc{Hd) = Pc to determine the representative pressure in the 
corona, and note that Pc y^Pd- Typically in the corona the represen- 
tative electron and ion temperatures (Tg^c and Jp c) are not equal. 
The corona is optically thin, and so we take the representative val- 
ues of pressure, density and temperature in the corona to be related 
by 



(9) 



Note the factor 2 in the denominator of eq. (|9]l. It ensures that for 
a one-temperature plasma (Jp c = T^^c) the correct value for the 
gas pressure is recovered (cf. eq. |2]i. We assume a pure hydrogen 
plasma for the corona (number density of electrons equals num- 
ber density of protons). An extension to non-hydrogen plasmas in- 
volves factors of a few percent in Given the other approxima- 
tions made here, not least the simplification of the one-zone model, 
we can safely neglect these factors for the purposes of this paper. 

We use the same viscosity prescription in the corona as for the 
disc (see eq.|4]i, adjusted to the corresponding values in the corona. 
The torque in the corona is given by 



Gc 



3 7 

-nVcpcHcR^c , 



and the heating rate in the corona is 

Qt = IvcPcHc^l 



(10) 



(11) 



The scaleheight of the corona may be comparable to the radial dis- 
tance from the black hole, i.e. Hc/R ~ 1. Then the coronal flow 
is pressure supported and deviates from Keplerian rotation. In the 
course of this paper however we keep the assumption of Keplerian 
rotation for the corona as well. As long as Hc/R does not greatly 
exceed unity, there are only factors of about two involved. They 
are however comparable to the likely errors we make in the vertical 
averaging and thus we neglect them. 

The heating in the corona is distributed onto the electrons and 
protons. The frac tion of h eat going into electrons and protons is not 
very well known. lOruzinovl ( 1 1 998h shows that if the magnetic fields 
are near to equipartition, high radiative efficiency is achieved and 
most of the viscous heat goes into the electrons, while for low mag- 
netic fields the heat goes mainly in the ions and low radi ative effi- 
ciency , as assumed by ADAF models can be reproduced. lOuataertI 



d 19981) . on the other hand, argues that there is damping at the pro- 
tons Larmor radius and hence most of the vi scous energy heats the 
proton s. This picture has been reconciled bv lOuataert & Gruzinovl 
(Il999l) . who show that for equipartition magnetic fields most of the 
energy heats the electrons, while for a small magnetic field the en- 
ergy primarily goes into the ions. They find a transition for = 5 
to Pb = 100, where (3b is the ratio of gas to magnetic pressure. The 
uncertainty then is introduced by uncertainties how Alfven waves 
are converted into whistlers on the protons Larmor radius. In the 
light of this discussion, we define our "fiducial" heating model for 
the corona by heating ions and electrons according to their partial 
pressure, P^^c and Pg.c, i-e- we take 



and 



2p,c Qc p 



Qlc = Qt^. 



(12) 



(13) 



respectively. This, together with the Coulomb collisions, intro- 
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duced below, leads to a preferred ion heating for a two-temperature 
plasma, and an equal share for a one-temperature plasma. 

Electrons and protons are coupled through Coulomb colli- 
sions. They transfer heat from the pro tons to the electrons at a ra te 
Qep. We take the rate as calculated bv lstepnev & GuilbertI (Il983h . 

The electrons cool radiatively by bremsstrahlung and other 
cooling processes (atomic line cooling etc.). While there are cool- 
ing rates for different metallicities and physical conditions available 
in the literature (e.g. . Suth erland & Dopita 1993) we limit ourselves 
here solely to electron-proton brems strahlung. We use the recently 
recalculated rate A^rems (7;,c, Tp.c) by'MaidEooS). The recalcula- 
tion accommodates the previously overlooked fact that for a two- 
temperature plasma the protons can contribute a significant fraction 
of the bremsstrahlung due to their plasma speed. For proton tem- 
peratures in excess of Jp c > "^p/"^e7e,c the proton thermal speed 
is larger than that of the electron and correspondingly most of the 
bremsstrahlung is produced by the deceleration of the proton rather 
than the electron. This calculation however strictly is only valid for 
the non-relativistic regime, although it could be generalised to in- 
clude this if required. 

Soft photons from the disc cool/heat the corona by Compton 
scattering at a rate (e.g. Rvbicki & Lightman 1979) 



Ae,c,C = Admax(Xes,Xes) 



max 40e,160^ 



(14) 



where Aj is the soft photon source (see eq. [3), Xgs = KgsPc^c 
the electron scattering optical depth (Kes = 0.3 cm^ g~^), 9e = 
^7e,c/('MeC^) and Teff^d the effective temperature of the disc. Both 
an Compton thick corona (Xgs > 1 or slightly relativistic electrons 
(0e > 1/4) lead to a strong increase in the cooling rate and therefore 
try to keep the electron temperature at Sg ~ 1 /4, i.e. Tg f« 1.5 x 10^ 
K, or lower. The first term in the square brackets correspond to 
Compton cooling, while the second reflects Compton heating. Note 
that the use of the effective temperature in (fT4l l for Compton heat- 
ing is valid only if the disc emission dominates. Then the inverse 
Compton temperature equals the effective temperature of the disc. 
In general, Tgff d should be replaced by a term proportional to th e 
flux weighted energy per photon (e.g. iRevnolds & Wilmsll2000h . 
For the cases presented in this paper, the electron temperature in 
the corona is always larger than 10^ K, i.e. Compton heating never 
plays an important role. 

For the corona we assume a local two-stream radiation field, 
i.e. half of the radiation (Bremsstrahlung and Compton radiation) 
is directed upwards and leaves the system, while the other half is 
directed downwards. If the disc exists, this part of the radiation field 
heats the disc at a rate 



2heat,d 



1 



(Ae,c,C + Abrems(7'e,c, Tp^c)) 



(15) 



where a is the albedo of the disc. For a = all radiation is absorbed 
while a = 1 corresponds to a perfect reflector. iHaardt & Maraschil 
(Il99lh estimate the albedo to be between 0.1 and 0.2 by Monte 
Carlo simulations. Without any significant loss of accuracy, we set 
the albedo a = 0, so that all the radiation impinging on the disc (if 
it exists) is assumed to be absorbed. 



2.3 The transition layer 

We now need to consider the structure within the transition layer 
between the disc and the corona. We first assume that the pressure 
across the transition zone is constant, i.e. we set 



(16) 



and hence in the transition layer of width AH (We essentially need 
and assume AH <^ (//(j,//c)) only density and temperature vary. 

In the transition layer heat is exchanged by thermal conduc- 
tion (see Sect.jBj. This process leads to mass and energy exchange 
between disc and corona by condensation and evaporation. The 
amount of mass flowing through the transition layer is assumed 
to be sufficiently slow that it does not affect the hydrostatic equi- 
librium and sufficiently fast that it can be treated as a stationary, 
i.e. time-independent, process, which only depends on the actual 
physical state of the corona and disc. 

We write for the ene rgy transport in the transition layer, fol- 
lowing Hegelman^&McKee (1990) 



d_ 



5k{Tp + Je 
' 2 2fAm„ 



+ ?cond = ^ 



(17) 



where rhz is the mass flow rate per (disc) surface area, ^cond is the 
conductive flux of heat given in ( IB4b and and q~ are volume 
heating and cooling rates as opposed to the vertically integrated 
heating and cooling rates used above. If > 0, then the disc evap- 
orates into the corona, and if riiz < the corona condenses into the 
disc. 

The characteristic length scale for eq. ( 117b is the so-called 
Field length (.Field. 19651 i.e. 



max (^+ , q~ 



(18) 



It can be shown that for a one-temperature plasma a stationary 
corona (accretion rate independent of radius), where the viscous 
heating dominates, leads to a ratio between the Field length and the 
height of the corona of 



67lKc, 



ond 



'T)T 



Ma^kT/ium^) 



(19) 



If we assume that the luminosity L of the disc corresponds to 
the accretion rate M via L = rjMc^, where c is the speed of light and 
r] = 1/ 12 the efficiency of converting rest mass into energy, we can 
scale the accretion rate and luminosity in terms of the Eddington 
luminosity and accretion rate Lpdd and Mndd- We get for the ratio 



^ = 1.7.10-3 



V^Edd/ 



109 K 



^ (20) 



Note that the result is independent of the mass. For R = 300 7?s and 
Te = Jp = 10*^ K and M = lO'^Mndd we get Xp/^c ~ 0.12. Hence 
the assumption is justified, at least with respect to the corona. 

The volume heating rate of the disc can be reduced to an ex- 
pression only containing pressure and rotation frequency by using 
the hydrostatic equilibrium (eq.[T]i and the dissipation rate ([6} to get 



-aPdQ.V2. 



(21) 



and similarly the dissipation rate of the corona (using eq.lTjandfTTTi 



Ql 

He 



-aPrD. 



He 



Hd + 



(22) 



Note that for a thick corona (He ^ H^) the volume heating rates 
show the same proportionalities (<=< aPD.). For a thin corona on top 
of a thick disc {He <^ H^), the volume heating rate is reduced by a 
factor {Hc/Hd)^~. 
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Figure 1. Schematic picture of our disc-corona sandwich, starting from the midplane (z=0) upwards. The disc has a half-thickness of H^, while the coronal 
layer has a thickness of Hq. The average pressure in the disc is approximated by the pressure in the midplane, and the pressure in the corona by the pressure 
at the base of it, Pc- For simplicity we plot the pressure as a linear function of the height. At the disc-corona interface the pressure at the base of the corona and 
the top of the disc is matched smoothly. We indicate the presumed temperature structure in the thin transition layer (TL) in the small inset. The pressure there 
is constant, but the temperature changes from the effective temperature of the disc Tgff ^ to the electron and proton temperature of the corona (rp ^ and Te,c), 
including a (possible) transition to a two-temperature plasma, if the ion temperature exceeds T* . Note that the length and pressure values are not to scale. 



The cooling terms can be calculated accordingly. For Comp- 
ton cooling we get the volume cooling rate from (fT4l i 



A, 



c,e,C 

He 



KespAd max 4ee,160 



(23) 



where we assumed a Compton thin transition layer. Given the geo- 
metrical thinness of the layer (see eq. I20t compared to the corona, 
this approach is justified. The volume cooling rate ( |23l l is strictly 

valid only for an entirely optically thin corona^ 

We now need to solve Ml\ . We follow Begelma n & McKe3 
([l990) and Rozanska & Czernv (2000b) and multiply by ^cond 
integrate over the transition layer to get 



m-, - — I ^cond az 



Hi dz \ 2 l/Limp 



1 



?cond (-^^d + ) ^ ?cond (-^^d ) ' 



is zero. If we now define the Field length using the result of 
iRozaiiska & Czemvl JioOObh . then we get 



(^TL ^Tl) ^cond'^Z 



(25) 



(24) 

The conductive flux at the boundaries of the transition layer 



Since the conductive flux contains a temperature gradient, the in- 
tegrals over the width of the transition layer in (l25T l can be trans- 
formed in an integral over the corresponding electron and proton 
temperatures. Note that this definition of the Field length differs 
from the one in dlSt . but for the relevant limiting case (coronal 
electron and proton temperature much larger than the disc values 
and either heating or cooling dominating the other) they are equiv- 
alent. For heating equals cooling the Field length formally goes to 
infinity, but then physically there is no mass flux due to conduction 

We are only interested in an average mass flux in the transition 
layer. Thus we take ritz = const, and hence can put ritz outside the 
integral on the right hand side of ( 124b . The knowledge of the Field 
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length ( 125) allows us to estimate an average conductive flux by 
setting dT /dz ~ T fkp where we assume that the coronal electron 
and proton temperature is much larger than the disc temperature. 
We finally find an estimate for the condensation/evaporation rate 
riiz to 

J'Ttif d (^TL ~ ^Tl) ^ond {^Te + ^epdTp) 2 IjUnip 

(26) 

where we set ^ep = m^/m^. For the heating we use the viscous 
volume heating rate ( 122b . while for the cooling we use the vol- 
ume rate for Compton cooling ((23) and Bremsstrahlung from 
^brems(7e,c5 7p,c)> which we divide by the corona scale-height to 
get the volume cooling rate, accordingly. 

The sign of mz crucially depends on the net effect of the heat- 
ing and cooling balance in the transition layer. If the cooling dom- 
inates the heating in the transition layer, then hot matter from the 
corona condensates into the disc (mz < 0). If the heating dominates, 
disk matter evaporates into the corona {m^ > 0). 

We show a schematic picture of the disc-corona sandwich in- 
cluding the transition layer temperature structure in Fig. [T] In the 
inset we show the assumed temperature structure of the transition 
layer. The proton temperature is assumed to rise monotonically 
from the photosphere, i.e. the effective temperature of the disc, to 
the coronal value. At some critical value T* however the proton 
electron coupling becomes very weak. Then electrons and protons 
start to have different temperatures. We set this critical temperature 
to T* =3-10^ K. For the evaluation of mz according to eq. ( |26t we 
need to take this consideration into account. For a one-temperature 
corona, we safely can neglect the integration over Jp, since it only 
enters with a factor of nie/nip. For a two-temperature corona, we 
need to integrate from Teff^d to Tp = Tg = T* and then until Tp c and 
Tg^c, respectively. 

We tabulate parts of the integral for selected values of Tp^c, 7e,c 
and Tgff d and use interpolations in these three variables to calculate 
thz in order to save computational time. 

2.4 Truncation 

At some point in the evolution of the disc-corona sandwich the disc 
or the corona can cease to exist at some radius. The determination 
of the exact transition to the corona or disc only state however is 
not straightforward. We discuss the two cases below. 

2.4.1 Disc truncation 

Here we take the simple approach that whenever the surface den- 
sity of the disc falls below a fraction of the surface density of the 
corona at that point, we assume the disc does not exist any more. 
We assume the disc ceases to exist if the surface density of the disc 
is smaller than 80 per cent of the coronal surface density. The disc 
is "reborn" if the disc surface density exceeds 90 per cent of the 
coronal surface density. 

The corona a priori never becomes opaque to absorption and 
does not reach values to become optically thick to electron scat- 
tering either. Hence our disc evaporation criterion approximately 
coincides with the point where the electron scattering optical depth 
of the disc, Xgs = Kes^c becomes lower than unity. Then our as- 
sumption of optical thick radiative transfer (eq. [3j cannot be used 
any more for similar reasons and different methods need to be used. 

If the electron scattering optical depth of the disc becomes 




Figure 2. Schematic picture of the truncated disc. The region between 
and Roui is filled with hot coronal gas only (a sample coronal column at 
R* is shown in pink), while the outer and inner parts still consist of an 
optically thick disc with a hot corona above and below. The coronal gas in 
the truncated region can only be Compton cooled by soft photons supplied 
by the outer and inner disc (shown in blue). 



lower than unity, then the high-energy radiation of the corona can- 
not be scattered or absorbed either and the heating dlSb cannot oper- 
ate any more. The disc becomes translucent. To which extent there 
can be an optically thin, but cool layer of gas, is not obvious. In 
the course of this paper we simply assume that it cannot exist and 
these parts of the accretion flow are completely filled with hot coro- 
nal gas. 

The truncation criterion used in this paper is somewhat unre- 
alistic, since even long before the disc becomes optically thin to 
electron scattering, it is optically thin to absorption and then it can- 
not thermalize internally produced or incoming radiation. The disc 
will heat up under such conditions and eventually become part of 
the corona. Our assumption however still can be used, since the ac- 
tual transition from a disc-corona sandwich to a corona only state 
takes place over a very small radial distance. Assuming a more re- 
laxed condition (i.e. truncate the disc if its surface density is smaller 
than 10^ times the coronal surface density), the results are still very 
similar. 

The hot corona-only phase in the middle has to be cooled by 
soft photons coming from the disc further out and/or, if it still ex- 
ists, from the inner disc as well. We show the situation schemati- 
cally in Fig. |2] The outer disc supplies a total rate of soft photons 
of 

5soft,out=47i/' ar^f d^JT? (27) 

For convenience we integrate to oo, i.e. the outer boundary of the 
disc, since oT^^^ oc and thus the integral only depends on the 
value of the effective temperature close to Rtr. At a radius R* < Ru, 
only a fraction of / = 27l7?*//(/?*)/(27i(/?tr -R*)^) passes through 
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the corona, i.e. 



f 



R*H{R*) 



(28) 



On the way from Ra to R* the amount of soft photons is reduced 
due to scattering at radii R with R* < R < 7?tr. Because the corona 
is optically thin to electron scattering in vertical direction, the same 
holds in radial direction (since He ~ R). Hence this dilution is not 
very significant and we neglect it here. S ince Xes ^ 1> the dilution 
factor, exp(— Xes) is close to unity. At R* the fraction / of the soft 
photon rate 5soft passes through the area 4nR*H{R*), i.e. we have 
the soft flux at R* 

^ "^SoftjOUt 



soft,out 



4%{Ri,-R*)^ 



(29) 



The electron scattering optical depth in the interval 

[R* -AR/2,R* +AR/2] can be approximated by 



Xes(^*)=KesPc(^*)A7?. 



(30) 



We calculate the Coulomb cooling rate at R* by multiplying by 
4^7c,e/('?JeC^). Division by AR leads to the corresponding volume 
cooling rate and multiplication with He to the final cooling rate 

>2 2 ^sP"^ 



i '-'soit,out "-^c,e ^ Tj 1 \ 

Ae,c,C,out = T^sPc^c (->!) 

Similar arguments hold for the Compton cooling which comes from 
the inner disc. Assume the inner discs ends at R^^ > 3R$ , then the 
ratio of the Compton cooling at R^^ < R* < /?tr is 



Ae,c,C,in _ f ~ R* \ 5'soft. 



A, 



e,c,C,out 



R*-Ru 



'^soft,o 



(32) 



The first factor on the right hand side changes from values much 
larger than unity for 7?* R^^ to values much smaller than unity for 
R* — > /?tr- The second term, the ratio between the soft photon rates, 
is proportional to Rm/Rtr (Here we assumed again Jgff d ■x R~^). 
Hence the contribution of the soft photon flux from the inner disc 
is, as expected, only large close to the inner disc. For the calculation 
we retain this contribution and write 



A, 



e,c,C,m 



^soft,in 



l^esPc^c 



(33) 



2.4.2 Corona truncated 

At some point in the evolution of the accretion flow, the cooling 
of the corona can exceed the heating and the scale-height of the 
decreases while the density increases. This is especially true for 
the parts of the corona further away from the centre. There cool- 
ing mechanisms in addition to those considered in this paper, like 
atomic cooling, make the corona thermally unstable and lead to 
runaway cooling. The cooling can be halted by Compton heating 
(cf. eq. [14] for Tgff d > r^g). Then the coronal matter effectively 
becomes part of the disc. It also may be halted by the onset of ra- 
diation pressure as the optical depth to absorption approaches or 
maybe even exceeds unity. The non-negligible optical depth will 
also cause photoionisation heating, since the soft photons of the 
disc get absorbed by the collapsing corona. The examination to 
which extent photoionisation can stabilise the corona is likely to 
be important in some circumstances but is beyond the scope of this 
paper. If at some radius the corona completely collapses and be- 
comes part of the disc, then these parts of the disc are still heated 
by the hard radiation coming from nearby disc annuli where the 
corona still exists. 



3 THE RELATIVE IMPORTANCE OF THE CORONA 
AND THE DISC FOR THE TOTAL LUMINOSITY 

To give an indication of the disc structure we are likely to be aim- 
ing at, we now estimate what fraction of the luminosity can be pro- 
duced by the different parts of a disc and corona sandwich. For 
simplicity, to make a rough estimate, we here neglect radial advec- 
tion and conductive losses/gains. The results of these estimates give 
a strong hint as to the radial extent of the corona and the disc, and 
we discuss some of the results from such considerations which are 
to be found in the literature. 

The stationary, i.e. time-independent version of eq. ( 143b is 

e;t+2heat,d=Ad , (34) 

and for the corona (combined for electrons and protons) 

2^ =Abrems+A;;^ ,. . (35) 

If the disc has an albedo a, then the corona heats the disc at a 



-a: 



■e,c,C 



The remainder of the 



rate 2heat,d = 1/2(1 - a) (^Abr, 
luminosity created in the corona escapes and equals the observable 
luminosity 



l+a 



Abrems A 



e,c,C 



(36) 



The luminosity of the disc is simply Ld = A^ . The coronal luminos- 
ity can be written as Lc = (1 +<3)/(l — fl)2heat,d- Hence the fraction 
Lc/Ld is 



U/LA = {l+a)/{\~a) 2heat,d /Aj 



(37) 



It is evident from eq. ( |34t . that 2heat,d is always smaller than A^ . 
Thus the maximum luminosity ratio for a passively heated disc 
can be Lc/L^ = (1 +a)/ (l — a) k, 1.35 where we set a = 0.15 
(iHaardt & MaraschJl99lh . Such a ratio applies for a passively only 
heated stationary disc. In all other cases the ratio will be lower. 

Coronal luminosities much larger than disc luminosities (rep- 
resentative for the low/hard state) can be produced in two obvious 
situations: either the disc is truncated and the inner accretion flow 
solely consists of hot, coronal gas, being cooled by soft photons 
from the outer disc, or we view the system along specific lines-of- 
sight and take into account relativistic effects. 

In the case of X-Ray binaries there is mo unting evidence that 
the disc is truncated in some states. [Esin et al.l ( 11997i) gives an uni- 
fied picture for black hole accretion flows. They use the accre- 
tion rate as indicator for the different states. For the low-hard and 
quiescent states they find a truncated disc, where the inner disc 
is filled with an ho t cor onal adv e ction dominated accretion flow. 
iHoman et all ( l200lb and iBellonil ( |2004|) introduce HID (hardness 
intensity diagram) as a phenomenological tool to distinguish dif- 
ferent states of X-Ray binaries in a unified scenario. Fender et al] 
(l2004 include contributions from a jet in this scheme. In their 
low/ hard state they assu me a truncated disc. 

iBarrio et"aD (l2003b discuss the geometry of the accretion flow 
in Cyg X-1. They fit the spectra in the energy range from 3-200 
keV with two models. One model utilises magnetic flares above an 
untruncated disc and the other assumes an truncated disc with an 
inner hot accretion flow. For the strong reflection in the magnetic 
flare case they do not find a break in the spectrum at high energies 
while with the weakly illuminated, truncated disc model they easily 
can fit the 3-200 keV spectra. 

iKubota & Donel ( Eto4) discuss the very high-state geometry 
of XTE J 1550-564. They calculate the luminosity Ld of the disc 
emission and show that in the very high state the luminosity is not 
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compatible with oc T^, where is the effective temperature of 
the disc. They conclude that the simplest explanation might be a 
smaller effective area where the disc emission is emitted, i.e. a trun- 
cate d disc. 

jPone & Gierlinskil (l2006h show that in XTE J1650-500 the 
broad iron line often attributed to extreme relativistic smearing can 
be well explained in the truncated disc geometry. Then the broad 
iron line is formed in an outflowing disc wind and the smearing is 
significantly reduced. 

For AGN, the situation does not look that dissimilar. However 
observable black-body components only can be expected from the 
very lower-mass AGN, since only then the effective temperature is 
accessible to present X-Ray satellites. 

Ptak_et aL (200 4) present observations of NGC 3398. They 
fail to find any signs of reflection from an optically thick 
disc. Thus they conclude that the disc is truncated at roughly 
100-300 Schwarzschild radii. Similar conclusions are found b 



and for the corona 



IStarling et all (l2005t) in the case of NGC 72 13. lzhang & Wul (Eopi 
use a truncated disc model to fit observations of the Seyfert galaxies 
NGC 5548 and NGC 4051. They find a disc truncated at 17. . .70 
and 700 Schwarzschild radii, respectively. Models using reflection 
(therefore assuming an optically thick disc all the way down to the 
last stable orbit ) however seem to fit the spectrum equally well (for 
NGC 4051, see lPontietall2006h . 



4 THE TIME-DEPENDENT EQUATIONS 

Having introduced the various physical considerations we require, 
we are now in a position to introduce the time-de pendent equations 
for our model. For the energy equations we followlMav er & Pringld 
( l2006h . We use the mass per unit area (i.e. the surface density Y) 
and the internal energy e per unit mass in the disc and corona as 
independent variables (for the corona we consider the internal en- 
ergy for the electrons, ee,c and protons, c, separately). All other 
variables can be determined by solving eqns. dAlb . (IA3b and (IA4b 
in Appendix lAl 



4.1 Continuity equations and angular momentum transport 



For the continuity equations we use 
aid 1 3Md 



dt ^ 2nR dR 



and 



dt 



1 aMc 



(38) 



(39) 



where and Zc are the surface density and and Mc the accre- 
tion rates in the disc and corona, respectively, is the mass trans- 
fer due to thermal conduction (see eq. [26l l. riiz > corresponds to 
disc evaporation, while < corresponds to the condensation of 
the corona into the disc. 

We calculate the accretion rates from the angular momen- 
tum transport equations where we assume an additional source/sink 
term for mass gain/loss. In combination with the continuity equa- 
tions ( 138b and ( |39b this source/sink term cancels out, since we 
take/put mass with the same lever arm a s angular momentum i s 
taken/put in (cf. the situation for R = Ra in lMaver & Pringlel2006l) . 
Thus we get for the disc 



Md 



d{R^aa) 



dR 



dR 



(40) 



Mr. 



d(R^ac 



dR 



dGc 
dR 



(41) 



where Hd and £lc are the respective rotation frequencies which we 
take to be Keplerian, i.e. Q.d = ^c = \J GM/R^. 

4.2 Conductive heat transport 

In Sect. |2.3| we calculate a mass transfer rate owing to the imbalance 
of heating and cooling in the transition layer (eq.|26b. 

If the hot corona condenses into the disc {rhx < 0), the residual 
heat is transferred to the corona at a rate Scond = — "^z^c- If the disc 
evaporates into the hot corona {thz > 0), the corona is heated at a 
rate ^cond = "^z^d- Thus 



5, 



cond 



— mz^c if mz < (corona condenses) , 
mz^d if «iz > (disc evaporates). 



(42) 



ed and ec are the specific internal energies per unit mass in the disc 
and corona (cf. eqns. lAll and lA2b . 



4.3 Energy equations 

The energy equation for the disc is 

^ (Id^d) + ^d^d = ^ (Md^d) — ^ — 

g^Vdd;-T- d d 2nRdR^ 2KRdR\pdJ (43) 

+ (2d + 2heat,d ~ ^d ) ~ "^cond , 

where Ld, ed, Pd and Hd are the surface density, specific internal en- 
ergy per unit mass dAlb . Pressure (|2]l and scale-height in the disc. 

and Aj are the local heating and cooling rates (cf . eq. [6] and 
[3j, 2heat,d the heating by the impinging radiation of the corona (cf. 
eq. [15} and S'cond the heat source/sink term for conduction ( 142b . 
The first two terms on the right hand side represent the radial ad- 
vection of energy and 'pdV'-work, while the others represent local 
source/sink terms. 

The energy equation for the corona is split into one for the 
electrons and one for the protons. We write for the protons 



dt 



Id,. . 
2kR dR ^ ' 

+ (e^c-eep-Abr, 



Pp,c d 



InR dR [ pc 



+ 5, 



cond" 



(44) 



and 



|(Ic.e 



I + Pe.cHc 



2kR dR^ 



' e,c " 

27iRdR 



rems.e 

•^cond " 



ee,c 

ec 



A 



e,c,C 



(45) 



where Jlc, ^p,c/^e,c> Pp,c/Pe,c and He are the surface density, specific 
internal energy per unit mass ( IA2b . Pressure ^ and scale-height of 
the corona. The pressure and internal energy is split into the respec- 
tive fraction for protons and electrons. Q^^ and Q^^ are the heating 
rates for protons and electrons, respectively (see eqns. llllll2l and 
[T3I 1. 2ep the electron-proton collision rate and Abrems,p and Abrems,e 
the electron-proton bremsstrahlung rate owing to the contributions 
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of protons and electrons, respectively. ^ is the Compton cool- 
ing/heating term (ea.ll4b. The first two terms on the right hand side 
of eq. ( 144b and ( 145 1 represent the radial advection of energy and the 
'pdV'-work, while the others are local source/sink terms. 

Note that 2ep is an interaction term and is treated as such in the 
energy equations. It conserves energy since the Coulomb collisions 
transfer energy from the protons to the electrons. The combined 
energy equation for the corona, the sum of eq. (l44l l and ( 145 1 does 
reflect this. The conductive energy transport conserves energy as 
well as is evident from the sum of all three energy equations. The 
same applies for the conductive mass exchange in the continuity 
equations. 



5 NUMERICAL SETUP 

Here we discuss the numerical scheme, the initial setup and the 
boundary conditions used in this paper. As we discussed above, the 
numerical results presented here are severely limited by the numer- 
ical resources required. We use a logarithmically equidistant grid 
with typically = 30 points starting from R = 3RstoR = 3000 Rs- 
All variables except the coronal and disc accretion rate Mc and 
are defined on the grid point. The accretion rate then is defined in 
between the grid points as is evident from eq. ( 140b and ( 141b . 

We take initial values for our simulation from a stationary 
model where we neglect radial advection and conduction. The (ini- 
tial) accretion rates in the corona and disc are fcM and (1 — fc)M, 
respectively. 

We integrate the system of time-dependent equations with a 
one-step Euler scheme and treat the advection terms in a first-order, 
upwind donor cell procedure. We have written our scheme so that 
we conserve mass and energy to machine accuracy. 

At the inner boundary we impose a zero-torque condition 
(Gd(^in) = Gc{Rm) = 0). This is equivalent to setting £c(^in) = 
^d(^in) = 0- All matter and energy flowing through this inner 
boundary is lost from the system. It is either added to the black 
hole or lost in a jet. No advection of energy and mass from within 
the inner radius is allowed. 

At the outer boundary we set the accretion rate M = for 
both the disc and corona. Just inside the outer boundary at the last 
disc cell, we supply mass to the optically thick disc at a fixed rate 
Mq = (1 —fc)M. Energy corresponding to the internal energy of the 
last disc cell is put into the last disc cell at a rate (1 — /c)Med- No 
mass nor energy is supplied to the corona. Thus Mq is the only pa- 
rameter characterising the system. The accretion flow then evolves 
self-consistently within the model and mass is transferred between 
the two phases accordingly. 

As one of the computationally more intensive parts of the sim- 
ulation is the calculation of the evaporation/condensation rate ( |26b . 
we use tables for essential parts of the integral and use interpola- 
tions in the course of the simulations. This enables us to calculate 
the evaporation/condensation rate with sufficient accuracy without 
spending too much computational time on the calculation. 

The computational time involved in the simulations is consid- 
erable. Since the geometries of the disc and corona are very dis- 
similar (the disc is geometrically thin while the corona is geometri- 
cally thick), the viscous timesc ales are vastly different. The viscous 
timescale is (e.g. |Pringlelll98lh 



which in units suitable for the presented models leads to a value of 



Xv = 1.4- 10"^ s- 



a \-i 



M 



0.1/ VlOMoy/ V^s / \R 



r^vr^i . (47, 



For the disc usually the ratio between scale-height and radial dis- 
tance from the black hole H/R ^ 10~^, while for the corona 
H /Rk. 1. Thus the ratio of the viscous timescale of the disc and 
corona at a given radius is about 10^. Thus it takes 10^ times longer 
for the disc to reach a (possibly) viscous equilibrium than for the 
corona. The thermal and hydrostatic timescales, however are simi- 
lar. 

Numeri cally the viscous ev olution allows for time steps 
smaller than dBath & Pringl3ll98lh 



A/ < - 



1 R^R 



2 llaQ-KH^ 



(48) 



For a = 0. 1, // f« 7? and AT? = 57? with 6 < 1, it is evident, that the 
maximum allowed time step for the viscous evolution is close to or 
even below the dynamical timescale. 

The dynamical timescale at the outer boundary at 7? = 3000 R^ 
is about 230 s, i.e. we need to integrate about 2 • 10^ times the dy- 
namical timescale of the innermost orbit to reach the dynamical 
timescale and another 10^ times to calculate for one viscous time 
at the outer boundary at /? = 3000 R^ . We then arrive at about 
2.3 • 10^ s to be calculated in real time, using about 2-10'^ time 
steps. 

Given the different timescales involved in this problem, an im- 
plicit solver would be ideally be suitable. We however chose to use 
an explicit solver since implicit solvers usually have problems with 
switching on/off processes (here mainly the mass transfer of eq. [26] 
which switches signs depending on the cooling and heating in the 
transition layer) given their large extrapolation capability. The dis- 
advantage of the explicit solver then is the large number of time 
steps needed to evolve the disc. 

In order to save computational time, we can make further ap- 
proximations. Formally all five time-dependent equations consti- 
tute a coupled system for the time-dependent variables ^d> 



Xv 



1 



aQ.K{H/RY 



(46) 



e^^c and ee,c (eqns.|38Jl39]|43]|44|and|45]l. All the other variables 
at a given time can be calculated considering hydrostatic equilib- 
ria (eq.[T]i and (|8j and the respective expressions for the pressure (|2]l 
and (|9]l and energy ( IAll and lA2l i. This system of 5 equations can be 
reduced to a system of 3 equations (see eq. I A 1 1 and I A3 1 - lA4l i but 
then it only can be solved numerically. However with a few simpli- 
fications a (computationally expensive) numerical solution can be 
avoided. It is evident that since generally <^ He, can safely 
be neglected in (jsjl. As we only consider only a disc with very low 
accretion rates, we also can neglect radiation for the disc pressure 
and internal energy (second term on the right hand sides of eqns.|2] 
and lAlb . By applying these approximations, we can analytically 
solve the five equations and hence avoid unnecessary numerical it- 
erations. The same applies to the terms involving H in both the 
coronal (eqns. I44landl45]l and disc ( 143b energy equation. A change 
dH can be expressed in terms of the surface density and internal 
energy of the disc and corona (see eqns. |A5| and |A6b . In the course 
of this paper, we choose to neglect this term in order to save compu- 
tational time. Physically this is justified since the corona and disc 
only influence each other through this term if the pressure in the 
disc and corona are close to each other or the scale-height of the 
disc and corona are similar (cf. eqns. IA5l and lA6b . Then however 
the disc or corona ceases to exist and other measures need to be 
taken (cf. Sect. 12.4b . In all cases other than that the coupling in- 
volving the H terms is negligible. 
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5.1 Disc develops a hole 

Our recipe for the truncated part of the disc are as follows. If the 
surface density of the disc falls below 80 per cent of the surface 
density of the corona, the disc does not exist any more. The corre- 
sponding cell however is allowed to advect mass and energy from 
the neighbouring still active disc cell. The disc is "reborn", if the 
"dead" disc surface density has refilled up to 90 per cent of the 
coronal surface density. "Dead" disc cells do not have thermal con- 
duction, produce no torque, are neither heated by the corona nor by 
viscous dissipation. The corona in this parts of the flow is cooled 
by soft photons from the outer and inner (as long as it exists) disc 
parts (cf. Sect. |2.4t . The innermost cell of the outer disc is defined 
as the cell where most of the mass is evaporated in the corona. 
The outermost cell of the inner remnant disc still loses mass to the 
next "dead cell" for angular momentum conservation. If the surface 
density, i.e. mass, of the disc in this cell is sufficiently high, then 
the disc is reactivated again. As the surface area of the disc annuli 
increases proportionally to R^, the outwards advected mass hardly 
ever manages to reactivate a disc cell. 

5.2 Corona fully condenses into disc 

This case does not occur in the model presented in this paper. If 
the corona did fully condense into the disc, then the corresponding 
torques and accretion rates in this parts of the corona decrease to- 
wards zero. If the scale-height of the corona shrinks below some 
limit (Physically we should take the point when the optical depth 
with respect to absorption becomes larger than unity), this part of 
the flow then is in a disc only state. Mass and energy influx from 
either side into the corona is still allowed. The corona starts to exist 
again if the optical absorption depth of the corona has shrunk below 
a certain threshold. 



6 RESULTS 

In the calculation we describe here, we consider an accretion disc 
model for a 10 M© black hole accreting at 10~^ ^Edd- We describe 
the time-evolution below. 

In Fig.[3]we show a schematic view of a typical sandwich disc 
in the initial stage of the time-evolution. We show relevant actual 
physical quantities 0.1, 20, 50 and 100 kiloseconds after the start of 
the simulation in Figs. 1417 1 

We start with a stationary disc and corona extending over the 
total computational domain. To compute the initial disc and corona 
structure we neglect radial advection, and we take the fraction of 
accretion occurring in the corona as /c =0.1. 

Initially (Fig.© in the outer regions the disc begins to evapo- 
rate into the corona, whereas in the inner regions the corona begins 
to condense into the disc (see Panel 1 in Fig. [3). Since the disc loses 
a lot of matter in the outer parts by evaporation into the corona, the 
corona temporarily starts to flow outwards (Fig. 14)1). Most of the 
matter is in the disc, and the corona is two-temperature inside a 
radius of 7? 100 7?s. 

Although the corona close to the black hole has got a very 
small surface density in the initial model, as the outer disc how- 
ever evaporates mass into the corona, this matter is then transported 
inwards through the corona. Further in, the coronal accretion rate 
decreases as part of the coronal matter condenses into the disc and 
fills up the corona to surface densities of about 0.3 . . . 1 g cm~'^. 
This process happens very quickly, i.e. on the viscous timescale 



of the corona. The increase in coronal mass in the inner regions 
leads to an increase in the coronal luminosity. The inner corona 
itself condenses into the disc and increases the disc surface den- 
sity (compare Fig.|4j5 and[5j)). Since most of the matter condenses 
very close to the actual transition radius between the evaporating 
disc and recondensing corona (which starts at around R k, 300 /?s 
and increases eventually to around R ^ 1000 R$), it takes about a 
viscous timescale of the disc at that radius to increase the surface 
density of the inner disc. This timescale is 10^ times longer com- 
pared to the viscous timescale of the corona. 

At the transition radius (where the flow changes from an evap- 
orating disc to and recondensing corona), the disc starts to flow out- 
wards (see Panel 1 in Fig.[3]and Fig.[5}i). The accretion flow in the 
disc is then divided into three zones. In the inner zone the disc is 
moving inwards, in the intermediate zone it is moving outwards, 
and all the accretion is carried by the corona. In the outer zone the 
disc is moving inwards. At the boundary of the intermediate and 
the outer zone, despite the fact that the disc flow is towards this 
radius, the surface density of the disc declines rapidly. This comes 
about because most of the mass flowing into this region from both 
sides is quickly evaporated into the corona. An empty ring, or gap, 
develops in the disc at this point. This gap consists of coronal gas 
only and is Compton cooled by soft photons from both the outer 
and inner disc remnants. 

Once this gap has formed (Fig. [6]l, the inner disc is accreted 
and the gap widens. A decline in the disc scale-height is already 
visible in Fig. [Sj;, but the truncation of the disc is fully evident in 
Fig. [6] The inner disc disperses. This state is shown schematically 
m Panel 2 of Fig. [3] 

Most of the inner disc gets accreted while a small part flows 
outwards (see Fig.[6ji). The luminosity of the disc declines. Eventu- 
ally we are left with an accretion flow where the inner part consists 
of a corona only and the mass supply for this corona is produced 
by the outer remnant disc (cf. Panel 3 of Fig.[3]and Fig.[7]i. 

We show the luminosity evolution of the corona and disc in 
Fig-ID Initially, the luminosity of the corona increases very quickly 
(caused as we described above by a filling of the inner parts of the 
corona on the short viscous timescale of the corona at the transition 
from the two- to a one-temperature plasma), but cannot exceed the 
disc luminosity as long as a disc exists (see Sect.[3]l. On the much 
longer, viscous timescale of the disc the disc luminosity increases 
(causes by a steady increase of surface density of the inner disc), 
before the luminosity drops with the development of the gap. Once 
the inner disc is completely accreted, the dominant part of the radi- 
ation comes from the hot corona. The disc becomes steady-state at 
this point, and no significant further evolution occurs. 



7 DISCUSSION 

The results of these computations indicate that a disc around a 10 
Mq black hole, accreting at M = 0.001 M^^d relaxes to a state 
in which the outer disc is a standard accretion disc with a small 
corona, but the inner regions (R < 1000 Rs) consist of only hot 
coronal gas which has evaporated from the disc. 

7.1 Thermal conduction and the transition from 
condensation to evaporation 

These calculations basically agree with the results of 
iRozatiska & Czemvl ( l2000bh who concluded that for a sta- 
tionary disc there would be net condensation of the corona in the 
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Figure 3. Schematic view of a typical sandwich disc in the course of the time-evolution. We only show one side of the sandwich. The sandwich is symmetric 
with respect to the horizontal axis. The aiTows indicate the mass flow within the corona (light red) and within the disc (blue) as well as the conductive mass 
flow between the two phases. The flow radially comprises of three zones: In Zone II (one-temperature corona) the disc evaporates leading to the formation of 
a corona. The evaporated mass flows inwards within the corona. In Zone I (two-temperature corona) the corona condenses into the disc. Both the disc and the 
corona accrete in this zone. At about the transition radius R,^ disc flows outwards (Panel 1). After some time, the disc is truncated and the hole is filled with 
coronal gas (Panel 2). The inner disc quickly gets accreted and the system is left with the situation as shown in Panel 3. 



inner disc and net evaporation further out. This comes about for 
the following reasons. 

The calculation of the evaporation/condensation rate (eq. [26l l 
shows that the sign of the term {q^i^ — q^i} determines conden- 
sation and evaporation, while the Kcond term given its strong tem- 
perature dependence determines the magnitude of the condensa- 
tion/evaporation. This strong temperature dependence also ensures 
that for the heating and cooling mechanisms considered in this pa- 
per the mass transfer rate mainly depends on the upper boundary of 
the integration, i.e. the coronal values for the temperature. 



To get an idea what causes the evaporation/condensation in the 
two zones shown in Fig. [3] we assume that the only heating is the 
aPc heating M2\ and Compton cooling ( |23l l and bremsstrahlung are 
the only cooling mechanisms. 

Within the corona the system is usually in thermodynamical 
equilibrium or quickly relaxes back to equilibriumm - faster closer 
to the black hole than further out. 

In the two-temperature plasma (zone I) Compton cooling 
dominates. If we calculate the difference between volume heating 
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Figure 4. Some physical quantities in the time-evolution of the accretion flow 100 s after the start of the simulation: Panel a shows the electron and proton 
temperature in the corona and the temperature of the disc. Panel b,c, and d show the surface density, the ratio (H/R) and the accretion rate for disc and corona, 
respectively. Note that we normalized the accretion rate to the rate we feed the disc in the outermost grid cell. The disc is evaporating in its outer parts very 
strongly and leads to an outflowing corona beyond a critical radius (see Panel d). 



and cooling, we get 

- %,c = Pc ttci^c - 8KesA/mp — f- , (49) 

\ + t^sC I 

where we have neglected Compton heating and the relativistic gen- 
erahsation in i23l . Within the corona, this term is about zero (as- 
suming thermal equiUbrium). To get the evaporation/condensation 
rate, we need to integrate the difference over the range of temper- 
atures from the disc to the coronal values. If the corona is a two- 
temperature plasma, then the fraction T^/ (Jg + Jp) is small, but the 
closer we get to the one-temperature plasma within the transition 
layer, the bigger it becomes. Then the difference ( |49l l becomes in- 
creasingly negative. The corona condensates onto the disc. 

In zone II we have got a one-temperature plasma and 
Bremsstrahlung dominates the cooling. We need to calculate (For 
simplicity we use the classical Bremsstrahlung result, i.e. qc brems °^ 

p2ri/2) 

it - Cbrems = ("0^0 " , (50) 

where C is an arbitrary constant of proportionality. While the dif- 
ference is zero in the corona, the further down we integrate towards 
the disc temperature, the smaller gets T and the larger the second 
term in the difference. Cooling dominates heating in the transition 



layer and thus the corona condenses into the disc. This is in contrast 
with our results presented in Section [6] The outer parts of the disc 
evaporate because the outward advection in the initial phases of the 
disc evolution lead to a imbalance between heating and cooling. 
Heating dominates there and hence the disc evaporates. 

7.2 Physical processes neglected in this model 

We now discuss some of the additional physical processes which 
we have neglected here, but which will need to be taken into ac- 
count if the model is to be developed further. 

7.2.1 Other cooling mechanisms 

Further out other cooling mechanisms than Compton cooling and 
bremsstrahlung will be more efficient and eventually lead to a con- 
densation of the corona. We neglect atomic cooling and only allow 
for bremsstrahlung and Compton cooling for both the corona and 
the transition layer. The temperature dependence of bremsstrahlung 
and Compton cooling ensures that the evaporation/condensation 
rate (cf. eq. [26} only depends on the upper limits of the integral 
involved. If we needed to include cooling processes other than 
bremsstrahlung (i.e. atomic cooling etc.), then the calculation of 
the evaporation/condensation rate would be much more difficult. In 
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Figure 5. The accretion flow 20 ks after the start: The corona production is at its maximum and so is the accretion rate in the corona. Most of the corona is 
produced at the transition from the two- to one-temperature plasma. As a result the disc flows outwards in some parts. Note that the ratio H /R for the disc 
starts to decrease. The labelling is the same as in Fig.|4] 



the transition layer the gas then can have a significant absorption 
optical depth which in turn leads to photoionisation processes and 
significant radiation pressure. Photoionisation and radiation pres- 
sure also influence the corona if it starts to collapse and may lead 
to an outflow (See Sect. |7.2.2] l. The increase in cooling may also be 
responsible for changing the sign of rh^ from evaporation to con- 
densation in the outer parts of the disc and thus will stop the transi- 
tion radius to move outwards. 

With the widening of the inner gap a lack of soft pho- 
tons develops and the electron temperature reaches temperatures 
in excess of 10^" K. This invalidates our use of non-relativistic 
Bremsstrahlung rates. However given the other underlying approx- 
imations in the model, the missing factors of a few should not influ- 
ence the results too much. In a more realistic model, pair production 
will set in (e.g. Svensson 1982) and help to cool the corona. Pair 
production increases the electron scattering optical depth and thus 
also leads to a rise in the Compton cooling rate. In fact, many, if 
not all observations of AGN hint at the existence of a population 
of hot electrons in a very narrow range of electron temperatures 
(about 100 keV, which corresponds to 10^ K (for NGC 4151, see 
iBeckmann et alj|2005t) . These authors deduce an optical depth for 
the corona of x 1.3. This is much higher than the values encoun- 
tered in the model presented in this paper. 

The study of these cases is very important but requires further 
work. 



7.2.2 Wind/Outflow 

We have neglected the possibility of a wind/outflow throughout 
these calculations. Physically given the large ratio of H/R for the 
corona, the cor ona is on ly marginally bound within the gravita- 
tional potential jBlandfor d & Begelman 1999). Thus it may well be 
that mass, energy and angular momentum is carried away in an out- 
flow, jet or wind. This might then eventually lead to a reduction of 
H/R until a stable state is found. This would then be an additional 
cooling mechanism for the corona. Observations of low-luminosity 
AGN and X-Ray binaries in their low/ hard state indicate t hat they 
are dominated by large-scale outflows jFender et al. 2004). These 
outflows can be driven magnetically (e.g. Maver & Pringlell2006l) 
or by radiation pressure as well (cf. Sect. 17.2. l] i. 



7.2.3 Corona production 

Another important question to be answered is the production of 
the corona. Our model presented in this paper shows that once 
we start with a tiny corona and consider thermal conduction, the 
corona is produced mainly from the outer optically thick disc. How 
the corona came into existence in the first place is not clear. As 
described in Sect. 12.41 we have run simulations for a one -phase 
black hole accretion disc but allowing for the transition to an op- 
ti cally thin flow. The simu lations are similar to the ones reported 
in lMaver & Pringl^ (l200d) . except that we now use a prescription 
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Figure 6. The accretion flow 50 ks after the start: The disc develops a gap and the accretion flow in the disc is truncated. The inner disc disperses, i.e. parts 
flow outwards, while most of the matter is accreted. The outer disc evaporates and produces the corona. The coronal accretion flow in the gap is stationary, i.e. 
the accretion rate is radially constant (see Panel d) given its much smaller viscous timescale. The labelling is the same as in Fig.|4] 



bv lArtemova et alj ( 1 19961) for the optically thick-thin transition. As 
reported before, we still find limit cycles for the same range of ac- 
cretion rates. During the high-M phase, where the inner accretion 
disc is depleted, this part of the accretion flow becomes optically 
thin and the temperature rises to about 10^ K close to the black 
hole. This might be a promising mechanism to produce a corona 
and a possible explanation for the so-called very high (VH) or steep 
pow er law (SPL) state in X-Ray binaries (Remillard & McClintockl 
I2OO6 ). This optically thin coronal part (being geometrically thick) 
then can flow out-/inwards on top of the accretion disc. 



7.2.4 Extra shear 

In the course of this paper, we neglect non-Keplerian rotation for 
the corona. We here discuss some possible effects of the inclusion 
of non-Keplerian rotation into the model. These could play an im- 
portant role in determining energy dissipation and the redistribution 
of angular momentum (i.e. accretion) in the disc and the corona. It 
will be important to investigate the outcome of including these ef- 
fects in future computations. 

Since the corona is geometrically thick, radial pressure gradi- 
ents lead to a deviation from pure Keplerian rotation. Corona and 
disc thus rotate at different speeds. At their boundary (the transition 
layer) torques arise which lead to angular momentum exchange be- 
tween the two phases. The work done by these torques leads to 



an additional heating of the corona, in particular the ions. Further 
angular momentum exchange arises from the mass transfer due to 
thermal conduction as then both layer rotate on different speeds. 
We give here a rough calculation of the torques and dissipation rates 
which might arise from these processes and the effect of the extra 
torque on the accretion flow. 

The rotation frequency 7^ ^^k = \/ GM/E? in the presence 
of a radial pressure gradient is given by 



^2 , J_9^ 



(51) 



where Pc is the coronal density, the pressure and R the radial 
distance from the black hole. We approximate the pressure gradient 
by dPc/dR = P/Rxp with Xp = dlogPc/dlogR. Since Pc ~ pcC^,c 
and Cs,c ^ ^cHc, where Cs,c is the sound speed and He is the scale- 
height of the corona, we get 



'l~iH,/Rfxp 
If AD. = — Dc, then for small Hc/R < I v/e can write 



(52) 



Xp 



(53) 



AO. _ 1 /He 
Q^~^2 \ r' 

Thus coronal and disc flow rotate at different speed. The differential 
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Figure 7. The accretion flow 100 ks after the start: The inner disc completely dissolved by accretion and the inner parts of the flow consist only of a hot corona. 
This part of the corona is stationary. Further out there still is a remnant accretion disc producing the corona. The labelling is the same as in Fig.|4] 
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Figure 8. Luminosity of the disc L(j and corona Lc. Initially the coronal luminosity increases quickly for 1 ks. Then the corona condenses into the inner 
disc leading to an increase in the disc luminosity (cf. Fig.|4]and|5). Starting at 40 ks, the disc luminosity declines as the gap widens. After about 70 ks the 
disc luminosity declines rapidly (The inner disc dissolves, cf. Fig.|6]and|7) and the corona dominates the luminosity. The small-amplitude variations in the 
luminosity arise from the coarse resolution and numerical fluctuations. The numbers annotated to the arrows on the time axis correspond to the figure number 
representing the respective state of the system. Note that we use a logarithmic scale for the time to differentiate the highlight the various timescales. 
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rotation leads to a free energy per unit surface area of about 

^Ec = \^^mR? ■ (54) 

We note that usually Id » and thus A£c ~ j^c (AOT?)^. This 
energy is released on a multiple of the local dynamical timescale, 
i.e. we set 



(55) 



where > 1 is a parameter and Q.~^ the dynamical timescale. 
The release of the free energy ( 154b leads to an additional heating of 
the corona at a rate 

Qcdiff ^^ = 1 {AORfa,/^' . (56) 

At ILd+l-c 

If we put (153b into the last equation, then we get 

Gcdiff = I i^KRf X^ne/£i . (57) 

The heating timescale for this process is 



'Cdiff,c 



2c,diff 



Sd 



He 



) Xp^fE.^;' ■ (58) 



Thus the heating timescale is longer than the thermal timescale Xc 
of the corona (Xc ~ (ttc^c)"^) for 



(59) 



The differential rotation leads to a torque which tries to speed 
up the corona and slow down the disc, i.e. the torque is 

Gdm-Qc4m/^c = l^^^{^KRf(^^) ^'p^^e! ■ (60) 

If only the torques given in eqns. dlOb and (|5]l would be at work 
in the corona and disc, respectively, then in a stationary Keplerian 
disc (neglecting mass exchange due to thermal conduction) matter 
flows at a speed 

1 fH\ 1 fH\^ 

"ft.visc = ;t«Cs ( ^ ) ~ t:^ 17 Xii • (61) 



R 



R 



As long as Q. decreases with increasing radial distance, i.e. = 
(3 log n)/ (3 log 7?) ^ —3/2 < 0, then matter always flows inwards 
as expected. If we include the effects of the extra torque ( [60t . then 
the contribution to the speed of matter is 



"R.diff '■ 



Gdiff 



(62) 



I7?a(2 + Xi2) 

For a near Keplerian rotating flow we approximate Xq, = —3/2 and 
put in eq. ( 160b to get 

2 / rr \ 4 



and 



1 Lr 



M/?,diff,d 



He 



Zd + \ R 



xI/e'^kR 



(64) 



for corona and disc, respectively. The relative contribution to the 
inflow caused by the viscous torques are 

^ = -^y^('^)( "^y^pfE.' (65) 

"R,visc,c 4ac Zd + Ic V / \R J 



and 

"R,diff,d 



1 



MR,visc,d 4ad Sd + r, 

Thus for 

1 Sd 



^\ fHc 



4ac Td + S. 
the corona starts to flow outwards and for 



fE,< 



1 



4ad Sd+Ic 



He 



xl 



Hc\ 



(67) 



(68) 



the inflow speed of the disc is dominated by the extra torque. 

For a typical disc and corona, we have %p — \, H^/R ^ \, 
H^/R K. Q.e K, I \/2 and Ed S> Tc, then we have 

• an outflowing corona for Je^ < 1 / (4cXc \/2) 

• and an inflowing disc dominated by the extra torque for /^^ < 
10V(4adV2). 

• Since we physically cannot resolve timescales smaller than 
the dynamical timescale, i.e. /^^ > 1, oCc has to be smaller than 
ttc < 1 /(4\/2) f« 0. 18 for both effects to work. 

If we include the effects of thermal conduction, then the limit 
on fE^ for the outflow of the corona decreases(increases) if the disc 
evaporates(corona condenses), since then the mass exchange trans- 
fers angular momentum as well. If the disc evaporates (m, > 0), the 
corona yields angular momentum at a rate of rhzR^AQ.. 

The extra torque can be produced by the poloidal component 
of the magnetic field, B^, in the disc which slightly overextends into 
the corona. The smaller speed transforms the poloidal field into an 
azimuthal field. It creates an twist leading to a azimuthal field of 
the order of 



Bit. 



Aa 



(69) 



after one orbit. Thus the torque produce by this process can be ex- 
pressed as 

B? AO. 



'mag 



R 



B^B, 
4k 



R- 



4% D.r 



(70) 



We furthermore can assume that the poloidal magnetic field 
strength B^ is essentially equal to the intrinsic disc magnetic field 
Bdisc- Thus we slightly underestimate B disc- If we further use 
^disc = 47iad^d dShakura & Sunvaevlll973h . then we can write for 
the torque in dTOl l 



'mag 



AD. 
^^ad/'d. 



(71) 



(72) 



For Q.C = Hk/ \/2 we get with eq. ( |53b 

Tmag = -V2(^Y^ ttdPdXP • 

Finally we can use Pd ~ PdC^ d ~ ^d^K-^d and arrive at 

-adV^(Y) (Y)'^'i(i^K^)'xp. (73) 



'mag 



The comparison of ( |72l i and (|6( 

1 Ic 

JEr 



allows to estimate /e^. We get 

2 / \ -1 



He 



8adV2 2:c + Id V R 
For typical values for the corona (Xp - 



Xp-if) (74) 
I, He ^ R) and typical 
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values for the disc (Zj k, lO^ ' ^^Zc and H^/R 10 ^) we get /^^ k, 
(10-2 . . . l)/(8V2ad). For ttd = 0.1 this leads to fE, ~ 10"^ . 1. 
Since we underestimate Sjiso another factor larger than unity ap- 
plies to the result in ( 174b . Hence our estimate of /^^ given the phys- 
ical picture of the magnetic torque discussed above is compatible 
with an outflowing corona. It however constrains /^^ to be not sig- 
nificantly larger than unity. 

Note that the inclusion of thermal conduction introduces a 
non-constant accretion rate Mc {R) and (7?) for the corona and 
disc, respectively, even for the stationary case. Then extra factors 
enter the estimations presented above. 

The physical process presented here might be responsible for 
driving a corona produced close to the black hole further outwards 
where it can be influence the spectrum and may condense into the 
disc (cf. Sect. lT23T l. 



8 CONCLUSIONS 

We present time-dependent simulations for a two-phase accretion 
flow around a black hole. We start with a flow which initially con- 
sists of a cool and geometrically thin but optically thick accretion 
disc, sandwiched by a hot and geometrically thick but optically thin 
corona. We show that as the disc evolves, the disc-corona sandwich 
vanishes inside a certain radius. A gap in the disc develops and the 
inner disc is fully accreted. In the final state, we are left with a 
coronal hot inner flow and a disc-corona sandwich further out. 

The general behaviour of the disc as described in this paper is 
not completely unexpected, since there is mounting evidence that 
the accretion flow around black holes at low accretion rates is trun- 
cated inside some radius (see Sect. [3]l. While we assume a stellar 
mass black hole of 10 solar masses, the result is likely also to apply 
to higher black hole masses, i.e. AGN. The result is in support of 
the unified scenario of lEsinetal.l ( ll997h . where the optically thick 
disc is truncated below a certain radius for low accretion rates and 
only reaches the last stable orbit for higher accretion rates close to 
the Eddington limit. 

For our model, we tried to keep the number of assumptions as 
small and the physics as simple as possible. But clearly there are 
still areas for improvement which we discuss in Section|7]and plan 
to include in the future work. The model presented here neverthe- 
less gives an indication that the standard disc at low accretion rates 
cannot exist close to the black hole. The physical reasons for this 
are discussed in Section lTTT] 
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APPENDIX A: TRANSFORMATIONS 

We define the specific energy per unit mass for the disc as 
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and the corresponding energy for the corona 
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Then, with the use of eqns. ([T]i, (|2]i, (|8]l, (|9]l and the boundary condition ( |16t we can solve for Pc and pd, the densities in the corona and disc 
and the temperature for given e^, ec,'LA = Pd^d and Ed = Pc^c- We arrive at a system of equations comprising of jAlb and 
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3 Pc VPd 2 pc/ 



Pd— + ^7d - TP'^^^ + o^^T = 0- (^4) 
/^WTp 3c \ 3 2 Pd 



The electron and proton temperature in the corona follow trivially from the specific energy of the electrons and protons (see eq. lA2b . The 
variation of and He in terms of the internal energies and Cc and surface densities Ed and Zc is given by 
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and 
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APPENDIX B: THERMAL CONDUCTION 

Thermal conduction occurs in media with a temperature gradient. Conduction induces an heat flow to cancel this temperature gradient. In a 
plasma, electron-proton collision occur typically after the particles travelled the Debye length. The Debye length is given by 

where Eq is the dielectric constant, e the elementary charge, m the reduced mass (for electron-proton systems it is essentially the electron 
mass), u the relative speed of electron and proton with respect to each other, n the number density of the plasma and InA f« 20 the Coulomb 
logarithm. 

For a thermal plasma, the average Debye length can be found by integrating over the corresponding isotropic Boltzmann distributions 
for electrons and protons, respectively. Thus 
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where Pg and Pp is the temperature of the electrons and protons and is the Boltzmann constant. Since both the electrons and protons have 
an isotropic velocity distribution, the effects of alignment and anti-alignment cancel and we do not need to treat any angle-dependent terms 
in the velocity u^. For further reference we define 

^.P = ^^. (B3) 
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Note that for a one-temperature plasma (Pg = Pp) the correction factor 1 + b,ep is close to unity and the classical Debye length is reproduced. 
The conductive flux of heat then can be described as the thermal average of nmu^Xi^idv /dx). We find 

^cond = (nmsU^B^'^ 
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where and Wp are the respective speeds of the electrons and protons. Together with the Debye length ( IB2t we find 

^cond y ^2 J 47ly/m^lnA y dx dx ) 

' V ' 

Note that we again find for a one-temperature plasma (Tg = Tp) the correction factor 1 +^ep is close to unity and moreover the electron and 
proton temperature gradient are the same and hence the proton temperature gradient contributions to the conductive flux can be neglected. 
The value of the parameter Kcond is 

Keond = 6.9-10-^re^ (l+^,p)^ergs-icm-iK-i . (B5) 

This is so mewhat larger than t he val ue quoted in IjBegelman & McKeeJ ( Il990l) . Their value (not containing the correction factor 1 + i^ep) 
is based on iDraine & Giulianil ( ITgsi) who use a higher value for the Coulomb logarithm and furthermore include effects of different ion 
species. Our value however is only about 20 per cent larger. Compared with the other approximations in the model presented in this paper 
we consider our value to be sufficiently accurate. 

The physical effect presented in this section is similar to Maver (2007). In both calculations the classical derivations are extended for 
a two-temperature plasma. Then the proton contributions cannot be neglected any more since their plasma speed can become comparable to 
the electron plasma speed. 



